THE X-RAY TRANSFORM FOR A GENERIC FAMILY OF CURVES 

AND WEIGHTS 

BELA FRIGYIK, PLAMEN STEFANOV, AND GUNTHER UHLMANN 

Abstract. We study the weighted integral transform on a compact manifold with boundary over a 
smooth family of curves F. We prove generic injectivity and a stability estimate under the condition 
that the conormal bundle of F covers T* M. 



Let M be a compact manifold with boundary. Let T be an open family of smooth (oriented) 
curves on M, with a fixed parametrization on each one of them, with endpoints on dM, such that 
for each (x,^) € TM \ 0, there is at most one curve 7^,5 € T through x in the direction of ^, and 
the dependence on (x, ^) is smooth, see next section. Define the weighted ray transform 



where w{x^£^) ^ is a smooth non- vanishing complex valued function on TM \ 0. We study the 
problem of the injectivity of /r,«; on functions on M. We impose no-conjugacy conditions on F 
that would guarantee that Ir,u; recovers singularities. Under that condition, we prove that /r,io 
is injective for generic L, w, including analytic ones, and that there is a stability estimate. This 
is a generalization of the X-ray transform arising in Computed Tomography which consists in 
integrating functions over lines provided with the standard Lebesgue measure. 

In |Mul| . Mukhometov showed that = 1, is injective, for 

any set F, provided that the curves 7 have unit speed, and is simple w.r.t. those curves. The 
latter means that for any two points x, y in (}, there is unique curve in 7 connecting them that 
depends smoothly on its endpoints. He later showed that this remains true if w is close enough 
to a constant in an explicit way. Stability estimates were also given. In dimension n > 3 there 
is no such known result for an arbitrary simple family of curves. On the other hand, if F is the 
family of the geodesies of a given (simple) Riemannian or Finsler metric, and w is close enough to 
a constant, injectivity and stability of Ir^w was established in |Mu21 lMu3l \AR\ \BG\ [R] . 

The transform /r,to is not always injective, even for simple F. An example by Boman [B] provides 
a smooth positive weight function w so that Ir,w fails to be injective in a ball in R^, where F consists 
of all straight lines. 

In the present work, we have incomplete data, i.e., we do not assume that we have a curve in F 
through any point in M in the direction of any vector (unless n = 2). On the other hand, we want 
{A^*7, 7 G F} to cover T*M, the latter considered as a conic set. We do not assume convexity of 
the boundary w.r.t. F. If F is a subset of geodesies of a certain metric, then some geodesies (not 
in F) are allowed to have conjugate points, or to be trapped, but we exclude them from F. On the 
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other hand, the result is generic uniqueness and stabihty, and Boman's result shows that this is the 
optimal one in this setting. 

Our approach differs from the works cited above and uses microlocal and analytic microlocal 
methods. Such methods are not new in integral geometry, see, e.g., [Gul IGuSll IGuS2l iGrUl 13 
BQ Q], but we use some recent ideas that led to new results in tensor tomography and boundary 



rigidity of compact Riemannian manifolds with boundary, see |SU31 ISU41 ISU5j . 



2. Statement of the main results 

Fix a compact manifold with boundary Mi such that M™* D M, where MJ"^* stands for the 
interior of Mi. We equip Mi with a real analytic atlas, where dM is smooth but not necessarily 
analytic. We will think of the curves 7 as extended outside M to Af™* so that their endpoints 
are m M|°*, and 7 n M remains unchanged. Different extensions will not change lT,wf as long as 
7 n M is the same. By 7^;^^, we will frequently denote the curve in T, if exists, so that x G 7^;^^, 
and 7^.^^ = /x^ at the point x with some /x > 0. We will freely shift the parameter on but not 
rescale it, so we may assume that x = 7^,^(0), then 7a;,^(0) = /x^. 

We want 7^.^^, for (x,^) G TM, to depend smoothly on (x,^), therefore in any coordinate chart, 
7 = 7^,^ solves 

(2) 7 = G(7,7), 

where = 7x,^(0) is smooth. The generator is only defined for |^| = |7x,^(0)| (in any 

fixed coordinates) but we can extend it for all In case of a Riemannian metric, for example, 
= — r^j(x)^^^', for l^lg = 1, and extended for all ^. The generator G determines a vector 
field G on TM that in local coordinates is given by 

(3) G = e'|, + O'(.,0|,. 

see also ([2^ . pUj) . The curves 7 G T are the projections of integral curves of G to the base, with 
appropriate initial conditions that reflect the choice of the parametrization. 

We assume that T is open with a natural smooth structure as follows. Fix any {"fit); < 
t < l^} G r, zizl± > 0, 7(/±) G M™* \ M, where we shifted the parameter t arbitrarily, and 
set xq := 7(0). Let H he a hypersurface in M™* \ M intersecting 7 transversally at xq, and let 
■^0 = 7(0). We assume that there exists a neighborhood U of (xq, Co) and a smooth positive function 
/i(x,C), (x,C) & U n H, with /i(a;o,i^o) = !> so that the integral curves of G with initial conditions 
(7(0), 7(0)) = {x,n{x,^)^), {x,0 e U n H, and interval of definition < t < 1+ belong to F 
(and in particular, the endpoints are in M{°* \ M). This makes F a smooth manifold; if H is given 
locally by = 0, then F is locally parametrized by {x',6) G R"""*^ x S^~^. We say that F is C^, 
respectively analytic, if G is G^, respectively analytic, on TMi, and for any such choice of G^, 
respectively analytic H, the functions fi are C^, respectively analytic, too. 

It is not hard to see that by duality, one can define Ir,wf for any distribution / G I5'(M™*) 
supported in M. 

Given x G M, we define the exponential map exp^{t, C), t G R, G TM\0, as exp^{t,(^) = "yx,^{t)- 
Note that exp^.(t,C) is a positively homogeneous function of order in the ^ variable, and in local 
coordinates, we can think that G S^~^. Then x = 7(0) and y = 7(^0) will be called conjugate 
along 7, if Dt^^ exp^(t, C) has rank less than n at (to, S,o), where = 7(0). It is easy to see that this 
definition is independent of a change of the parametrization along the curves in F (that we keep 
fixed). We would like to note here that (in a fixed coordinate system), the map v = t^ 1-^ exp^{t,£,), 
where |C| = 1, t G R, may not be C°°. In case of magnetic systems, for example, it is only G^ 
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while exp2.(i,^) is a smooth function of all variables, see [DPSUj . This requires some modifications 
in the analysis of the normal operator ([6]) , see section 14. 1[ 

It is clear that one cannot hope to recover any / from /r.to/; if there is a point in M so that no 
7 G r goes through it. We impose a microlocal condition that requires something more than that, 
we want any (x, Q) € T*M \ to be "seen" by some simple 7 G F. 

Definition 1. We say that T satisfying the assumptions above is a regular family of curves, if for 
any {x,Q € T*M, there exists 7 G F through x normal to C, without conjugate points. 
We call any 7 as above a simple curve. 

If F is not regular, one can give the following example of a non-injective Ir,w Let M be a 
subdomain with boundary of the sphere S*"^^ with its natural metric. Assume that M™* contains 
a pair of antipodal points a and b. Then any function that is supported in two symmetric to each 
other small enough neighborhoods A ^ a., B 3 b., and odd with respect to the antipodal map, 
integrates to over any geodesic in M. Not only Ir,wf with w = 1 does not determine /, it does 
not determine the singularities, either. For example, if / = Sa — ^b^ where 5a^b ^-re delta distributions 
centered at a and b, respectively; then /r,i/ = 0. 

On the other hand, one can see that Ir.wf, known for a regular family of curves resolves the 
singularities of /. Using analytic microlocal arguments, we also show that one can recover the 
analytic singularities, as well, if F is analytic. This allows us to prove the following. 

Theorem 1. Let F be an analytic regular family of curves in Mi and let w be analytic and non- 
vanishing in M . Then Ir,wf = for f G I?'(M™*) supported in M implies / = 0. In particular, 
It,w is infective on L^{M). 

To formulate a stability result, we will fix a parametrization of F. Let H he a finite collection 
of hypersurfaces {Hm} in M™* that are allowed to intersect each other. Then H may not be a 
hypersurface but is still a manifold if we think of each Hm as belonging to a different copy of M. 
Let H be an open conic subset of {{z,6) G TMi; z € H, 9 ^ T^H}, and let ±l^{z,9) > be two 
continuous functions. Let F(7Y) be the subset of curves of F originating from 7Y, i.e., 

(4) T{n) = {7.,e(i); i-{z, e)<t< i+{z, e), {z, e)en}. 

We also assume that each 7 G F('H) is a simple curve. 

We will fix a parametrization of a subset of F that is still regular. 

Given Ti as above, we consider an open set ?i' ^ TC, and let T{Ti.') (s T{Tl) be the associated set 
of curves defined as in (jj]), with the same l^. The restriction 7 G T(TC') C F(7Y) can be modeled 
by introducing a weight function a in TY, such that a = 1 on TC', and a = otherwise. It is more 
convenient to allow a to be smooth but still supported in Ti. 

We consider Ir,w,a = olr^wi or more precisely, 

(5) Ir,w,af = a{z,e) w{jz^e,iz,e)f{lz,e)dt, {z,9)en. 

Jl-(z,0)) 

Next, we set 

(6) Nr It 

Here the adjoint is taken w.r.t. a fixed positive smooth measure dS on 7Y; more precisely, we assume 
that in any local coordinate chart, dS := a{z, 9) dSz d9 on H, where dSz is the surface measure 
on H in the so fixed coordinate system, d9 is the surface measure on S^~^, and C°° 3 a > 0. 
Notice that dE is not invariant under a different choice of Tl and a coordinate system on it. On the 
other hand, injectivity of Nr^w,a is equivalent to injectivity of Ir,w,a, and the latter is equivalent to 
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injectivity of Ir^w restricted to suppa, see |SU3j . and this property is independent of the choice of 
fi and the coordinates on it as long as they parametrize the same set of curves. 

Theorem 2. 

(a) Let Ti.' <£: 7i be as above with r(H') C r(7-^) regular, and {G, fi, a,w) fixed. Fix a S C°° with 
TC' C suppa C Ti.. If It, w, a is injective, where T = T{H), then we have 

(7) \\f\\LHM)/C < \\Nr,u,,af\\HHAH) < C||/||l2(m). 

(b) Let H' (^H, a = be as above related to some fixed {Gq, fio,ao,wo). Assume that Iro,wo,aO 
is injective, where Tq = Tq{'H). Then estimate ^ remains true for (G, fi,a,w,a) belonging to a 
small neighborhood of {Gq, fio,ao,wo,a^), with a uniform constant C > 0. 

Remark In fact we need only regularity for w, a. 

We notice that above refers to different spaces. More precisely, fi, are considered in G'^iTi), 
while G, w are considered in C'^{TM). To define correctly G'^iTM), we fix any finite atlas on M, 
see also the remark in section [H 

Example (simple systems). Let M C be diffeomorphic to a ball, and let G{x,^) be a 
smooth generator on TM \ = M x R" \ 0. Fix a coordinate system on M. We can assume that 
G is defined on SM ^ M x 5""^ and extend homogeneous of order to all ^ 0. Set 

d-SM = {z£ dM; 6 ■ u < 0} 

where z^(z) is the exterior unit normal to dM. Then we define T as the set of all curves 7 = ^^,6 
that solve 

(8) ^ = G(7,7), 7(0) = 2, 7(0) = A(z,0)e, {z,0)Gd-SM, 

where A > is a given smooth function on MxS'^^^ with X[z, 6) = \'^z,e{^)\- Let 7^^^ be the maximal 
curves with those initial conditions. Assume that for any x G M, the map exp^, : exp~^(M) ^ M is 
a diffeomorphism depending smoothly on x. Note that this implies that all those curves are of finite 
length; for any x, y in M, there is unique 7 G F that passes through them, smoothly depending on 
X, y, and the curves in F have no conjugate points. As above, 7's are allowed to be directed curves; 
if X £ M™*, 9 € S"'~^ then the curves 'jxfi and 7a;, -6» are not necessarily the same. We also assume 
that Ml 3) M (meaning that Mj°* D M = M) is another domain diffeomorphic to a ball so that 
{G, A) extends smoothly there and satisfies the same assumptions. 
For a simple system as above, define 

(9) lG,x,n.f{z,0) = [ u;(7,,e(t),7.,e(t))/(7.,e(i))dt, {z,e) G d^SM,. 



One could also study subsets of curves as above. Let a be any positive C function on d^SMi, 
and set dS = a{z,9)\i'{z) ■ 6\dSzd9. Then 

Ig,x,^ : L\M) ^ L\d^SMi,dE) 

is a bounded map, and Ng,x,w = Iq \w^G,\,w is a well defined operator on L'^{M) that can be 
extended as an operator from L'^{M) to H^(Mi). Note that the factor \i^{z) • ^| in dS can be 
omitted since dMi is convex and M stays at a positive distance from dMi. If Mi = M, and if dM 
is strictly convex w.r.t. F, then that factor is needed to preserve the mapping properties of Nc^x^w, 
see |SU3j for the Riemannian case. 
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3. Injectivity OF /r,u, for analytic systems 

In this section we prove Theorem [TJ We denote by WFa(/) the analytic wave front set of /. 

Proposition 1. Let 70 G T 6e a simple curve. Let /r,u)/(7) = for some f € T>'{Mi) with 
supp/ C M and all 7 € neigh(7o). Let T and w be analytic near 70. Then 

(10) iV*7onWFA(/) = 0. 

Proof. We will choose first a coordinate system (x' , 2;") near 70 so that the latter is given by x' = 0, 
= t, t ^ [l~,l^] with some itZ^ > 0, and moreover, replacing x' = by x' = z, where z is a 
constant vector with \z\ -C 1, one still gets a curve in T (parametrized by t again, i.e., a unit speed 
line segment in the so fixed coordinate system). 

Fix a point po £ 70, and shift the parametrization of 9 t 1— s- 7o(t) so that po = 7o(0). 

Assume that po ^ M and that the part of 70 corresponding to ^~ < t < is outside M, too. 
Set X = expp^(t,9), where |^| = 1, t > 0, where the norm is in any fixed coordinate system near 
Pq. Then {t,9) are local coordinates near any point on 70 fl M because the 70 is simple. Since 70 
may self-intersect, they may not be global ones. On the other hand, there can be finitely many 
intersections only, and one can assume that each time 70 intersect itself, it happens on a different 
copy of M X R. More precisely, (t, 9) 1— > (t, expp^{t, 9)) C M x R is a codimension one submanifold 
of M for 9 close to ^0 = 7o(0) and t E (0,/^) by the simplicity assumption, and we think of any 
function / : M — > C as defined on that manifold. Therefore, without loss of generality, we may 
assume that 70 does not self-intersect. 

Write x' = 9', x^ = t. Then x are the coordinates we were looking for in 

U = {x; |x'| < e,r <t <l+] C Mi 

with < e ^ 1. They are analytic, since F is analytic. 

Fix xo € 7o, and G T*Mi conormal to 70. We need to prove that 

(11) (xo,e°)^WFA(/). 

By shifting the x"^ coordinate, we can always assume that xq = 0. Note that '■= 7o(0) = e^. 
Here and below, Cj stand for the vectors d/dx^ , and e^ stand for the covectors dxK 
Assume first that / is continuous in M and vanishes outside M. 

The arguments that follow are close to those in |SU5j . Set first Z = {x" = 0; |x'| < 7e/8}, 
and denote the x' variable on Z by z' . We will work with the curves t 1— > 7(z',o),(6",i)(^) defined 
on l~ < t < l~^, the same interval on which 70 is defined. Each such curve is in F for \9'\ <C 1 
because the latter is open. They all have endpoints in M™* \ M, and in fact, we modified a bit the 
endpoints of the interval of definition to make them constant (l^). We can do this, when e <C 1, 
and this does not affect integrals of / over them. 

Let Xn{z'), N = 1,2, . . . , he a sequence of smooth cut-off functions equal to 1 for \z'\ < 3e/4, 
supported in Z, and satisfying the estimates 

(12) [a"x7v| < (CA^)'"', ia|<A, 
see [Eil Lemma 1.1]. Set 9 = {9', 1), \9'\ < 1, and multiply 

Ir,wf (7{2',o),e) = 

by XNiz')e^^^ , where A > 0, ^' is in a complex neighborhood of {^^Y, and integrate w.r.t. z' to 
get 

(13) / / e'^''<'xNiz> (7{.',o),9(i), 7{.',o),e(i)) / (7(.',o),e(t)) (t) dt dz' = 0. 
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For \0'\ <^ 1, {z',t) £ Z X (/^,/+) are local coordinates near 70 given by x = 7(^',o),e(0- Indeed, 
if 9' = 0, we have x = {z',t). Therefore, for 9' fixed and small enough, {t,z') are analytic local 
coordinates, depending analytically on 9'. In particular, x = {z' + t9',t) + O(|0'|). Performing a 
change of variables in ([13]) , we get 



(14) / e^^"'(^'^')-«'a7v(x, 9')f{x) dx = 



for \9'\ <^ 1, VA, V^', where, for \9'\ <^ 1, the function {x,9') is analytic, independent of 

N, and non-zero for x in a neighborhood of 70, satisfies (112^ everywhere, vanishes for x ^ U; and 
0^(0,6*') = w{0,9). 

Without loss of generality we can assume that 

e = e^-\ 

Here and below, ej stand for the vectors d/dx^ , and stand for the covectors dx^ . 
We choose the following vector 9{£,) analytically depending on ^ near ^ = 

(15) 9(0 = fCl, . . . , j" + ---+^n-2 + U ^ ^ 

If n = 2, this reduces to 6'(^) = {—^2/^1, !)• Clearly, 

(16) 9(0 -^ = 0, ^"(0 = 1, ^(e°) = e„. 
Differentiate (fT5]) to get 

(17) = . = l,...,n-2, ^(e°)=0, ^(e°) = -e._,. 

In particular, the differential of the map S*"^^ 9 ^ 1-^ ^'(C) is invertible at ^ = = e^^^ . 

Replace 9 = {9\ 1) in (|14p by 0(^) (the requirement \9'\ ^ 1 is fulfilled for ^ close enough to 
to get 

(18) / e'^'^(-'^)a^(x,0/(3;)dx = 0. 



where (p is analytic in f7, and cln has the properties of ajsj above for ^ close enough to . In 
particular, 

ajv(0,0 =«^(0,^(e)). 

The phase function is given by 

(19) (^(x,e) = ^'(x,^'(e))-e'. 

To verify that is a non-degenerate phase in neigh(0, ^'^), i.e., that det (^^,^(0, ^'^) / 0, note first 
that z' = x' when = 0, therefore, {dz' / dx'){Q,9{^)) = Id. On the other hand, linearizing near 
x" = 0, we easily get {dz' / dx'^){id,9{i)) = -9'{i). Therefore, 

(20) <^,(o,e) = (e',-e'(6-0 = e 

by (fT6]l . So we get (^^^(O, ^) = Id, which proves the non-degeneracy claim above. In particular, 
X 1-^ (/3g(x,i^) is a local diffeomorphism in neigh(O) for ^ G neigh(^''), and therefore injective. We 
need however a semiglobal version of this along 70 as in the lemma below. 

Lemma 1. There exists (5 > such that 
for X £ U, \y\ < 6, |^ — ^'^l < S, complex. 
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Proof. We will prove the lemma first for y = 0, ^ = x' = 0. Since S,) = 0, we need to prove 
that the only solution to f^{{0, x'^),(^^) = in the interval /~ < < is x" = 0. 

We start with the observation that '/'(7o,(e'(5),i)(0) = 0- Differentiate the latter w.r.t. ^ at 
^ = ^0, t = x", to get 



where Ju{t) = d"fQ^Q{t) / dOy aX 9 = Cr, 



,n 



we get by ([I7D, (recaU that ^" = e""^), 



- 1, are "Jacobi" vector fields. Since (/?(x,^° 



(21) 



1, 

n 
n, 



where J" is the (n — l)-th component of J,^. Now, assuming that the l.h.s. of (j2ip vanishes 
for some fixed x" = to, we get that JJ!'^^{to) = 0, = 1, . . . ,n 



1. On the other hand, S := 
span(Ji(to), • • • , Jn-i(io)) is a hyperplane transversal to by the simplicity assumption. Therefore, 
for the unit normal to S, we have f„ 7^ 0. Hence, and e„_i are linearly independent, and the 
intersection of S and e:^_i is of codimension 2, and Ji(to)) • • • ) Jn-i{to) all belong there. Therefore, 
Jjy(to), = 1, . . . ,n — 1, form a linearly dependent system of vectors. The latter contradicts the 
simplicity assumption. 

The same proof applies if x' 7^ by shifting the x' coordinates. 

Let now y, ^ and x be as in the Lemma. The lemma is clearly true for x in the ball -6(0, ei) = 
{|x| < El}, where ei <C 1, because (p{0, is non-degenerate. On the other hand, V5^(x, ^) ^ ^^{y,C) 
forx€f/\S(0,ei),2/ = 0, e = C°- Hence, we still have (/^^(x,^) 7^ ^^{y,0 ^ small perturbation 
of y and ^. □ 

We will apply the complex stationary phase method [Sjj, see also |KSUt Section 6]. For x, y as 
in LemmaH and \ri - < 6/C, C » 1, 5 < 1, multiply (HH) by 

where x is the characteristic function the complex ball B(0,5), and integrate w.r.t. to get 



(22) 



3^^*(2''^'^'«)6jv(x,e,??)/(x)dxdC = 0, 



where is another amplitude, analytic, independent of N, and elliptic near 70 x {^^}, satisfying 
(fT2D. and 



-V9(y,0 + V5(2;,e) + 2(^-^)^- 



We study the critical points of 1-^ If y = x, there is a unique (real) critical point = and 
it satisfies '^^^s^ > at = ^c- For y ^ x, there is no real critical point by Lemma [H On the 
other hand, again by Lemma [H there is no (complex) critical point if |x — ?/| > 6/Ci with some 
Ci > 0, and there is a unique complex critical point if 1^; — y| < S/C2, with some C2 > Ci, still 
non-degenerate if (5 ^ 1. For any Cq > 0, if we integrate in (p2|) for |x — y| > S/Cq, and use the 
fact that |<I>^| has a positive lower bound (for ^ real), we get 



(23) 



\x-y\>5/Co 
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Estimate ([23|) is obtained by integrating N times by parts, using the identity 

as well as using the estimate (112p . and the fact that on the boundary of integration in ^, the e'^* is 
exponentially small. Choose Co ^ C2. Note that 9$ > for ^ G 5(suppx(- — and j] as above, 
as long as (7 » 1, and by choosing Cq ^ 1, we can make sure that is as close to rj, as we want. 
To estimate ([22]) for \x — y\ < S/Cq, set 

i{j{x,y,r]) := ^'l^^^^- 

Note that = — i(y — x) + 7/ + 0(6), and ipix, y,r]) = r/ • (x — y) + — + 0{6). The stationary 
complex phase method ISjl, see Theorem 2.8 there and the remark after it, together with ()23p . gives 



(24) / e^^'^(^'°)/(2;)S(a;,a;A)d2; = 0(iVe~^/^ + A"/2(C3iV/A)^), V7V, 

J\x-y\<S/Co 

where a = (y, rj), and S is a classical elliptic analytic symbol independent of A^. Moreover, the 
principal symbol ap{B){0,0,r]) equals w{0,6{r])) times an elliptic factor, and is therefore elliptic 
itself. Recah that w{0,6{(^)) = w(0,e„) / 0. Take so that N < X/iCse) < iV + 1 to conclude 
that the r.h.s. of ^ is 0(e~^/^). 
At y = X we have 

(25) ipy{x, X, rj) = -ipx{x, rj), ip^^x, x, r]) = ip^{x, rj), ip{x, x, rf) = 0. 
We also get that 

(26) 9V(2/,x,r?) > Ix-ylVC^, 

that can be obtained by writing y = x + h, and expanding tp in terms of powers of h up to 0{h^). 
Define the transform 

a I — > (5 = (a^., Va^v9(a)) , 
where, following a = {ax-,oi^)- This is equivalent to setting a = {y,r]), P = {yX)j where 
C = ^yiy^ f])- Note that C = ^ + 0{5), and at y = 0, we have C, = ri,hy (pO|) . It is a diffeomorphism 
from a neighborhood of (0,.^'') to its image, leaving (0, ^'^) fixed. Denote the inverse map by a{(5). 
Note that this map and its inverse preserve the first (n-dimensional) component and change only 
the second one. Plug a = a{j3) in (p^ to get 

(27) / e^^^(^'^)5(x,/3;A)/(x)dx = 0(e-^/^), 

J \x-ax\<&/Co 

for j3 S neigh(0, ^''), where iIj, B are (different) functions having the same properties as above, 
except that now tp satisfies 

(28) Vj/(3;,x,C) = -C, V'x(x,x,C) = C, V'(3;,x,C) = 0. 

By (SjI Definition 6.1], (l27|), ([28]), together with the eUipticity of B imply that 

(0,e°) 0WFa(/). 

Note that in [SjJ, it is required that / must be replaced by / in ([27]) . If / is complex- valued, 
we could use the fact that I(JR.f){'y) = 0, and /(5/)(7) = for 7 near 70 and then work with 
real- valued /'s only. 

If / is a distribution, then one can see that (jl4p still remains true with the integral in the x 
variable understood in distribution sense. The rest of the proof remains the same, except that 
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the cutoffs w.r.t. x fiave to be replaced by smootli ones. Tlie cliaracterization of WFa(/) in [Sj] 
Definition 6.1] is formulated for distributions, too. 

Tliis concludes the proof of Proposition [TJ □ 

Proof of Theorem [IJ The proof of Theorem [1] now follows immediately. By Proposition [U / is 

analytic in Mi and has compact support there. Therefore, / = 0. □ 

4. The smooth parametrix 



Under coordinate changes x x', G preserves its form, i.e., 

d I- , , d 



(29) G = e4j^+G'\x\a ^ 



and the transformation law is 

(30) G'Hx\i') = G^(x\^i'^)^ + 

V ox'i J ox'^ dx^ox'^ dx^ ax* 

This shows that the assumption G € C'^ is independent of the choice of the coordinate chart, 
and choosing a different finite atlas will preserve inequalities of the kind \\G — G\\q2^xm) ^ C'e by 
changing C only. 

We construct below a parametrix for Nr^w,a assuming that G, A, w, a are smooth. 

Proposition 2. NY^w,a an elliptic classical ^DO of order —1 in M"^*. As a consequence, there 
exists a classical pseudodifferential operator Q in M™* of order 1 so that 

QNr,^,af = f + Kf 

for any f G 2?'(M™*) with supp/ C M, and an operator K with a C^(M{°* x M{°*) Schwartz 
kernel. 

As a first step towards the proof of Proposition [21 we derive a formula for /p^^. Notice 
that the map TL x {1^,1'^) 3 {z,9,t) ix,v) G R" x S"'~^ given by x = exp^{t,6), v = 
dtexp^{t,9)/\dtexp^{t,6)\ is a local diffeomorphism. Indeed, fix {zo,9o,tQ), and let {xo,vq) be 
the corresponding (x, v). To find the inverse of that map, we need to solve 

ex.p^{—t,v) = z, —dtexp^{—t,v) =9, z £ H, 

for {z, 9, t) near {zq, 6*0, to)) so that {z, 9, t) = (zq, ^o, to) for (x, v) = (xo, vq). This can be done, since 
H is not tangent to any 9 such that (z, 9) G TL. Let J''(x, v) = d(z, 9, t)/d(x, v) be the corresponding 
Jacobian (depending on the choice of the local chart near x). 

Let (p e C^{n), f e C°°(M), and let wi € C^(rM{'^*) have smaU enough support that can fit 
in a coordinate chart that we fix. Then 



{Ir,n.„Md^ =111 a{z,9)wi{j,4t),%^g{t))f{^,,e{t))Hz,9)dtdS,d9 
a\x, v)wi{x, v)f{x)(j)^{x, v)J^{x, v) dx dv, 



where a'^{x,v) = a{z{x,v),9{x,v)), i.e., equals a, extended as constant along the curves ^z,eit)', 
and the meaning of (p^ is the same. Therefore, 

Ir,wi,a'Pix) = / a^ix,v)wi{x,v)(j)^{x,v)j\x,v)dv. 

J\v\=l 
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Let wq € C^(TM™*) be another function with smah enough support. Then 

Ir,wi,jT,wo,af{x) = {Jwi){x,v){a^wo){eyip^{t,v),dteyip^{t,v)) 

(31) ^5"-! J 

X f[ex.p^{t,v))j^{x,v)dtdv. 

The simpUcity assumption impUes that for any to / 0, and {xo,vo) belonging to the support of 
the integrand above, the map {t,v) i— » y = exp^{t,v) is a diffeomorphism from a neighborhood of 
{to,vo) to its image, and this is true for x in some neighborhood of xq. On the other hand, near 
to = 0, the map (t,u) ^ y = exp^(t, f) has Jacobian vanishing at t = 0, and the "true exponential 
map" = tv ^ y = exp^(t, v) is only a diffeomorphism, in general. By a compactness argument, 
given e > 0, one can cover M with finitely many charts, so that when x belongs to either one of 
them, one can split the integration in (j3ip into finitely many open sets that cover R\(— e, e) x 
In each of those integrals, perform the change of variables (t, v) ^ y = exp2.(t, v). Then we get that 
the l.h.s. of (jSip is an operator with a smooth kernel. The only contribution to the singularities of 
the kernel may therefore come from t G (— e,e). 

To analyze the contribution to (f3T]) from t G (— e,e), we proceed as follows, see also [DPSU] . 
The function m{t,v;x) = (exp^(t,t;) — x)/t is smooth, therefore 

(32) ex.p^{t,v) — X = tm{t,v; x), m{0,v;x) = X{x,v)v. 
We introduce the new variables (r, w) G R x S*""^ by 

(33) r = t|m(t, f ; x)|, uj = m(t,v;x)/\m(t,v;x)\. 

Then (r, lo) are polar coordinates for y — x = ruj in which we allow r to be negative. Clearly, (r, iv) 
are smooth at least for e small enough. Consider the Jacobian of this change of variables 

(34) J(x,t,T;):=det§^, 

d{t,v) 

computed with the same choice of local coordinates on 5"^^ for v and u) (and independent of that 
choice). It is not hard to see that J\t=o = A 7^ 0, therefore the map R x S^~^ B {t,v) 1— (r, w) G 
R X S^~^ is a local diffeomorphism from neigh(O) x S""^ to its image. We can decrease e if needed 
to ensure that it is a (global) diffeomorphism on its domain because then it is clearly injective. We 
denote the inverse functions by t = t{x,r,uj), v = v{x,r,ij). Note that in the (r, cj) variables 

(35) t = r/X + 0{\r\), v = oj + 0{\r\), %^,{t) = Xoj + 0{\r\). 
Another representation of the new coordinates can be given by 

r = sign(t) \e^p^{t,v) - x\ , to = sign(t)r — — r, 

|exp^(t, V) — x\ 

and 

(t, v) = exp~^(x + ruj) 

with the additional condition that r and t have the same sign (or are both zero). 

We return to (13ip . The paragraph after it shows that one can multiply the integrand by a smooth 
function x{t) so that X = 1 near t = and suppx is small enough; and the error is a smoothing 
operator. Then one can write, modulo a smoothing operator applied to /: 

Ir,wi,Jr,wo,af{x) = / x{t)Bix,t,v)f{exp^{t,v))dtdv 

(36) i 

xit)J^^(x,t,v)B(x,t,v)f(x + ruj)\. ^ , ^drda;. 
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where, see ([3T|) . 

(37) B{x,t,v) = {a^'Wi){x,v){a'^'Wo)(^exp^{t,v),dtexp^{t,v)^J^{x,v). 

4.1. Certain class of integral operators with singular kernels. Let U C R" be open and 
bounded. The integral representation (j36|) shows that we need to study integral operators with 
singular Schwartz kernels (with integrable singularity at x = y) of the class below, see also |DPSUt 
Appendix D]. 

Lemma 2. Let A : Co{U) — > C{U) be the operator 

(38) Af{x)= / / A{x,r,Lu)f{x + ruj)drdLO, 

with A G C°°{U X R X S^^^). Then A is a classical ^DO of order —1 with full symbol 

oo 

a{x,^) ~ ^afc(x,^), 

A:=0 

where 

afc(x,0 = 27r^- / d^A(x,0,uj)6^''Huj-ndio. 

Proof. Notice first that if A is an odd function of (r, cj), then Af = 0. Therefore, we can replace 
A above by A(,^ea{f,ijj) = {A{r,uj) + A[—r, —uj))/2. Next, it is easy to check that we can integrate 
over r > only and double the result. Therefore, 

oo 



(39) Af{x) = 2 / ^even(x,r,w)/(x + ru;)drdw. 

Consider now r, oj as polar coordinates for z = ruj, and make also the change of variables y = x + z 
to get 

(40) Af{x) = 2 [ A,,Jx,\y-x\,ji^)^^^^dy. 

J V \y - x\J \y - x]"^^^ 

Let 

7V-1 

(41) Aevenix,r,u;) = ^ Ae^en,k{x,uj)r^ + r'^RN{x,r,u) 

k=0 

be a finite Taylor expansion of Aeven in ^ near r = with > 0. It follows easily that 
2Aeven,A;(a^) ^) = Ak{x, uj)+{—l)^Ak{x, —uj), where klAk = d^\r=oA, and in particular, Aeven,k{x, oj)r^ 
is even w.r.t. (r, w). The remainder term contributes to ([^0]) an operator that maps L1^^^{U) into 
H^^^°{U) with some fixed Nq. To study the contribution of the other terms, write 

(42) Aven,fc/(X) = 2 j Ae,en,fc [x, ^3^) 1^ " dy. 

The kernel of -Agven.fc is therefore a function of x and z = y — x, with a polynomial singularity at 
y — X = 0, and it is therefore a formal ^'DO with symbol that can be obtained by taking Fourier 
transform in the z variable. Motivated by this, apply the Plancherel theorem to the integral above 
to get 



Aven,fc/(x) = {27t)-^ I e" «afc(x, 0/(6 d^. 
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where 



oo 

e-''''^<Ak{x,u;ydrdu 



(43) = 27ri*^/ Afc(x,a;)(5('=)(w -Odo;. 



In the third hne, we used the fact that Ae^cn,k{x,u>)r^ is even. Note that a^^x,^) is homogeneous 
in ^ of order —k — 1 and smooth away from ^ = but a distribution (in S') near zero. To deal with 
this, choose x ^ C'o" supported in |^| < 1 and equal to 1 near C = 0. Write a(x,^) = x(C)'^(^)0 + 
(1 — x(0)^(^) 0- The second term is a classical amplitude, while the first one contributes the term 

(44) Aven,fc(X*/) 

to (I42p that is smooth, as can be easily seen by making the change of variables z = y — a; in (j42p . □ 

Proof of Proposition\^ For x in a small enough neighborhood of a fixed xq, using a partition of 
unity {xj}, we can express Ny. yj^ci &s a finite sum of operators of the kind (|3ip. namely A^r,ui,a — 
w cJv w a with Wi = Xi'W- By the analysis following ([3T]) . their Schwartz kernels are smooth 
if we integrate outside any interval containing t = 0, and the only non-smooth contribution may 
come from terms of the kind ()36p . where wq and wi are replaced by some Wi and tfj. By Lemma O 
(j3ip is a classical \I'DO of order —1. Its principal symbol is given by 

0^0(2^10 =27r / A{x,0,uj)6{ll! ■ ^) duj. 



In case of (j36|) . A(x, 0, a;) is given by = J ^(x, 0, cj) (|a^tu|^XiXi) {x, w). Then ^-^^ is elliptic 
because w 0, and because given (x,^), there is w _L ^ so that a'^{x,Lo) ^ 0, and there exists i so 
that Xii^:^) > 0; and all other terms are non-negative. Therefore, A^r,«;,a is an elliptic ^'DO of 
order —1 in Mj°*, and the proposition follows. □ 

The next proposition is a standard consequence of the ellipticity of NY^w,a- See [SU3i. Theorem 2] 
for a similar statement in tensor tomography. In contrast to |SU3| however, we do not lose a 
derivative. 

Proposition 3. Under the conditions of Theorem\^ without assuming that Ir,w,a is injective, 
(a) one has the a priori estimate 



|L2(jvf) < C\\Nr^w,af\\m(Mi) + Cs\\f\\H~^(Mi): Vs; 

(b) Ker/r,M;,o is finite dimensional and included in C^{M). 

Proof. Part (a) follows directly from Proposition [2j Next, if / € Ker/r,«;,a, then (Id + K)f = 0, 
and K is a. compact operator on L'^{M), with smooth kernel. This proves (b). □ 
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5. Reducing the smoothness requirements 

In this section, we will reduce the smoothness requirements on F and the weight w, and will 
prove Theorem [2j 

We start with the observation that assuming that Ir,w,a is injective on L'^(M), then Nr^w,a '■ 
L'^{M) H^{Mi) is injective, also, see |SU3] . Then we get by Proposition [3l|b) and [TaH Propo- 
sition V.3.1] that 

(45) \\f\\LHM)<C\\Nr,^,e.f\\HHM,)- 

The second inequality in ([7]) is obvious. This proves part (a) of Theorem [2j 

In the rest of this section, we will perturb P, w, a and show that this will result in a small 
constant times ||/||l2(^) that can be absorbed by the l.h.s. above. We think of P as determined by 
(G, fi, a). Since Nr,w,a is a ^'DO that depends on P, w,a in a continuous way, if the latter belongs 
to C'^, A; ^ 2, the statement of Theorem [2]^b) follows immediately from what we already proved 
if there is replaced by C^, k ^ 2, see also [S\J3\ \S\J4\ ISU5] . Our goal here is to reduce that 
smoothness requirement. 

Proposition 4. Assume that G, fi,a,w,a are fixed and belong to C^. Let {G,jl,d,w,a) be 0{6) 
close to (G, fi, a, w, a) in C^. Then there exists a constant C > that depends on an a priori bound 
on the norm of {G, iJ,,a,w,a), so that 

(46) \\{Nr,^,^-Nr,^,c.)f\\HHM,)<Cm\\LHM) 

Proof. Assume now that we have two systems {G, fl,a,w,a) and {G, fi,a,w,a), as in the proposi- 
tion. Let Co be a bound on the norm of the first system. All constants below will depend on 
Cq. Let 5 be as in the proposition. To estimate the difference of those quantities related to the two 
systems, we will need the following comparison inequality for ODEs of Gronwall type. 

Lemma 3. Let x, x solve the ODE systems 

x' = F{t,x), x' = F(t,x), 

where F, F are continuous functions from [0,T] x U to a Banach space B, where U d B is open. 
Let F be Lipschitz w.r.t. x with a Lipschitz constant k > 0. Assume that 

\\F{t,x) - F{t,x)\\ <d, \ft£ [0,r], Vx G U. 

Assume that x{t), x(t) stay in U for < t < T . Then for < t < T , 

(47) \\x{t) - x{t)\\ < e^'\\x{0) - x(0)|| + ^ (e'^* - l) . 

Por a proof see |CLj . Note that the lemma can be used to compare the derivatives of x and x 
w.r.t. the initial conditions by differentiating w.r.t. the initial conditions first, and then applying 
the lemma. Since the curves 7 solve the equation ^ considered in the phase space, see ([3]), we get 
under the assumptions of Proposition HI 

(48) \\jx,v - %,v\\c^ + \\'yx,v - %,v\\c^ ^ ^"(5, 

where the norm is w.r.t. {x,v,t). The inclusion of t can be easily deduced from equation ([2]). 

Assume that G, /u, a, w, a are fixed and belong to C^. We will determine first the smoothness of 
the functions r{x,t,v) and u!{x,t,v) defined in ([33]) . Since 

(49) m(t,v;x) = / 7^,^,(s^)ds, 

Jo 
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we get that m and rh are functions of their arguments. By (j33p . we get that c^r and d^u), 
j = 0,1, are C^, also. In particular, the inverse functions t{x,r,uj), v{x,r,uj) are C^, too. On the 
other hand, J and in (|36p . (j37p are C^. Moreover, the difference of those functions for the two 
systems is 0{6) in the corresponding norms. 

Let us analyze first (|3ip in the case when the kernel there is multiplied by 1 — compare 
with (|36p . As explained in the paragraph following (j3ip . we perform the change of variables 
{t,v) ^ y = ex.p^{t,v) that is in our case, and its Jacobian is C^. Moreover, the Jacobians for 
the two systems differ by 0{6) in the norm. Then we get an integral operator with a kernel, 
vanishing near the diagonal x = y. Clearly, such an operator maps LF'^M) into C^(Mi). Moreover, 
the difference of two such operators, related to {G, fl,a,w,a) and {G, fi,a,w,a), respectively, has 
a norm 0{6). 

The more interesting case is what happens near the diagonal. To analyze this, we expand B in 
dSZD and J, see ^ as 

(50) B{x,t,v) = Bo{x,v) + tBi{x,t,v), J~^{x,t,v) = Jo{x,v) + tJi{x,t,v). 

The explicit expressions for Bq, Bi, Jq, Ji are listed below: 
Bo{x,v) = {{a'^fwiWo){x,v)j\x,v), 

Bi{x,t,v) = / bi{x, st,v)ds, where 
Jo 
d 

bi{x, t, v) = —{a^wi){x, v){a^wo) ( ex.p^{t, v), {dt exp^)(t, v)) j\x, v), 
Jq{x,v) = J^^{x,0,v) = \^^{x,v), 

1 



Ji{x,t,v) = J —^{x,st,v)ds. 

Notice that Bq,Bi, Jq, Ji € C^. Moreover, Bq,Bi, Jq, Ji differ by them by 0{6) in the norm. 
Then for A in ()38p . we get 

^^^^ A{x,r,uj) = x{t{x,r,uj)J~'^{x,t{x,r,uj),v{x,r,uj))B{x,t{x,r,uj),v{x,r,uj)) 

= : ^0(2;,^) + rAi{x,r,Lo). 

Here Aq(x,uj) and Ai(x,r,uj) are functions of all variables, and we have used the fact that 
t{x,r,io)/r is G^, too. As above, we get 

(52) \\Ao{x,uj) - Ao{x,uj)\\ci + \\Ai{x,r,uj) - Ai{x,r,Lo)\\ci <C5. 

Let ^1 be as in Lemma [2] related to Ai and rAi, respectively. Then the Schwartz kernel of 
see pUj) . is 2^o,even(a^5 '^)^"~"''''^) where we use the notation 

r = \x — y\, uj = [y — x)/r. 

Therefore 2AQ^evea{x,u:)r~'^~^^ has singularity of the type while the kernel of Ai has singu- 

larity of the type r"""*"^. To estimate the norm of A, we need to analyze the operator with 
kernel dxA. We get that formally, dxAo is an operator with a non-integrable singularity of the type 
r~", while 9a;^i,evcn IS an operator with kernel that still has an integrable singularity. Let now ^0,1 
be related to {G, p,,a,w,a). The contribution of Ai — A\ to ([Ml) is easy to estimate using (f52|) . 
The remaining question is whether 

(53) <C5||/||i2(^), fGL^iU), 
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where U is as in Lemma [21 and in our case, is a small enough open set in a fixed coordinate 
chart of Mi. We showed above that an operator with a weakly singular kernel, and dxAo is 
formally an operator with singular kernel. The continuity properties of the latter class are well 
studied, see e.g., [Stl IMP| . and the integration is understood in principle value sense. By the 
C alder on- Zygmund Theorem, see, e.g., |MP1 Theorem XI. 3.1], [St], a singular operator with kernel 
K{x,y) = Q{x,Lo)r~"' is bounded on L'^{U), if Q has a mean value in the u variable, and belongs 
to L°°{Ux', L'^{S^~'^)). Then the norm of that operator is bounded by C||0||, where the latter norm 
is in L°°{Ux; 

In our case, we start with an operator with weakly singular kernel 2^o,even 

(x,tj)r-'^+i that IS even 

w.r.t. uj, since ^o.evon it is independent of r. Therefore the x-derivative, if we differentiate the occur- 
rence of X in r and uj only, is an odd function of to. This makes the kernel dxi^AQ^Q^enix, uj)r~'^~^^) 
a singular odd one, up to a weakly singular kernel. Now \MP\ Theorem XI. 11.1] says that this is 
actually the kernel of dxAo, and by the Calderon-Zygmund Theorem, its norm is bounded 

by C||Ao||ci. We apply now those arguments to Aq — Aq with the aid of ([52]) . This yields ([53]) and 
completes the proof of the proposition. □ 

Proof of Theorem\M We already proved part (a) in ([^5]) . Combine estimate ([^5]) and Proposition [J] 
to get 

l|/|lL2(Af) < C\\Nr,w,af\\m{Mi) 

This immediately implies Theorem [2]^b). □ 
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